Eliminating the mean-field shift in multicomponent Bose-Einstein condensates 
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We demonstrate that the nonlinear mean-field shift in a multi-component Bose-Einstein condensate 
may be eliminated by controlling the two-body interaction coefficients. This modification is achieved 
by, e.g., suitably engineering the environment of the condensate. We consider as an example the 
case of a two-component condensate in a tightly confining atom waveguide. Modification of the 
atom-atom interactions is then achieved by varying independently the transverse wave function of 
the two components. Eliminating the density dependent phase shift in a high-density atomic beam 
has important applications in atom interferometry and precision measurement. 
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s-wave scattering, the potential is of the form V = 
(4-7T h na /m) 5 (r 12), where m is the particle mass, n the 
atomic number density, ri2 the relative position of the 
atoms, and a the s-wave scattering length. This leads 
to a density-dependent phase shift in the evolution of 
the condensate wave function. In the language of nonlin- 
ear optics this process is known as self-phase modulation. 
Since the exact density of the condensate is usually not 
perfectly known (especially in cases where it has been di- 
vided by a 'beam-splitter') this is a source of uncertainty 
that limits the accuracy with which precision measure- 
ments can be done. 

There is at first sight nothing that can be done to elim- 
inate this shift short of reducing the condensate density 
to a point where it becomes negligible, or of taking ad- 
vantage of Feshbach resonances, in which case three-body 
collisions appear to become a serious problem Q. The 
situation is quite different, however, for multi- component 
condensates. In addition to a self-phase modulation pro- 
portional to its own density, each condensate component 
experiences in that case an additional cross-phase mod- 
ulation, i.e. a phase shift proportional to the density of 
the other component p7|-pl|. As we will demonstrate, 
it is possible to engineer the environment of the BEC 
so that the phase shifts associated with self- and cross- 
phase modulation cancel each other, thus eliminating the 
density-dependent mean-field shift from the condensate 
evolution. As the property which governs the evolution 
of the condensate phase is the chemical potential, the 
elimination of the mean-field shift is equivalent to having 
a chemical potential which is independent of the number 
of atoms in the BEC. In addition, we show that it is pos- 
sible for at least one branch of the quasiparticle spectrum 
to be density-independent as well. Measurements which 
excite only the 'collisionlcss' branch are therefore insen- 
sitive to nonlinear phase shifts and hence of considerable 
interest in atom interferometry. 

To illustrate how this works we consider the simplest 
possible situation that can lead to the required effect, 
namely a two-component condensate described by the 
many-body Hamiltonian 



The promise of Bose-Einstein condensates |yj2|] and 
atom lasers HJlas sources for atom interferometry Q] 
based sensors [|3| results from their high brightness and 
coherence which leads to an increase in the signal-to- 
noise ratio as compared to conventional atom optics. Ad- 
ditionally, further enhancement of sensitivity might be 
achieved by taking advantage of nonlinear effects which 
occur in quantum-degenerate atomic fields. For exam- 
ple, feedback between optical and/or matter-wave fields 
can result in nonlinear instabilities [0,[lO|. Such positive 
feedback between optical and atomic intensity gratings 
has already led to the design of matter-wave amplifiers 
pl],p p^Jl^ | . The utility of such instabilities for nonlinear 
interferometry lies in the associated increase in sensitiv- 
ity that can be achieved by operating just above or just 
below the instability threshold. This allows a small per- 
turbation to produce a large change in the properties of 
the system. 

The transition to high density atomic samples carries 
a price, however, in that atom-atom collisions may in- 
troduce unwanted nonlinear (density dependent) phase- 
front distortions which limit sensitivity. While collisions 
are often viewed as a source of decoherence in quan- 
tum optics, in the regime of coherent atomic matter 
waves atom-atom interactions lead to coherent nonlin- 
ear wave-mixing and can therefore be manipulated using 
techniques inspired from nonlinear optics. In this Let- 
ter we discuss how the density-dependent phase shift can 
be eliminated in a trapped multi-component condensate. 
With applications in atom interferometry in mind we 
consider specifically the case of a narrow atomic waveg- 
uide such as those recently microfabricated on glass chips 
Jl4],ft5| , as these devices hold great promise for the devel- 
opment of integrated atom-intcrfcrometric devices. 

At temperature T ~ a single-component Bose- 
Einstein condensate is characterized by a scalar order 
parameter whose evolution is governed by a nonlinear 
Schrodinger equation (NLSE). At densities low enough 
that three-body collisions can be neglected the NLSE 
contains a cubic nonlinearity whose form is determined 
by the two-body collision potential. In the limit of 
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Here, Uj includes both the external trapping potential 
and the internal atomic energy for component j, and the 
constants g a , gb, and g x give the strengths of the non- 
linearities due to atom-atom collisions. In the language 
of nonlinear optics, g a and gb determine the self-phase 
modulation of the two components, whereas g x governs 
cross-phase modulation between them. 

The specific physical system that we have in mind con- 
sists of a condensate in a one-dimensional atomic waveg- 
uide with tight transverse confinement in the x-y plane, 
i.e., Uj(r) — Uj(x,y), while propagation along the z- 
dimension is free. The two components are two internal 
states (e.g., the Zeeman sublevels) of the same atomic 
species. We assume that they may convert into each 
other through a weak linear coupling. The ground state 
wave functions ipj (r) are hence found by minimizing the 
Hartree energy functional while holding the total num- 
ber of atoms fixed (although the individual atom num- 
bers may vary). This linear coupling is either sufficiently 
weak, or is turned off adiabatically after the steady state 
has been established, that we do not explicitly include it 
in the Hamiltonian. 

As ideally one would like to have a 'single-mode' wave 
guide, we assume that the transverse confinement is suf- 
ficiently strong that all atoms are 'frozen' in the ground 
state (pj(x,y) of the transverse potential Uj(x,y). This 
leads us to introduce the atomic annihilation operators 
(f>j(z,t) for atoms in the transverse ground state as 



4>j(z,t)= / dxdyip*Ax,y)^j{Y,t). 



(2) 



From the Hamiltonian ([!]) the Heisenberg equation of mo- 
tion for 4>j(z) is found to be 
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V^){z)h{z) + Vj\{z)4> k {z) fa(z), (3) 



where k ^ j and 

Vj = gj J dxdy \ip } (x, y)f ; j = a,b 
V x =g x I dxdy\ip a (x,y)\ 2 \(pb(x,y)\ 2 . 



(4) 
(5) 



We choose the energy reference such that E a = —5/2 
and Eb = 5/2. From these expressions it is immediately 
apparent that the effective phase modulation constants 



along the z-dimension can be modified by appropriately 
varying the transverse trapping potential Uj(x,y). 

Assuming that any perturbation present is too weak 
to excite the transverse degrees of freedom, the ground 
state wave function and quasiparticle spectrum of low- 
lying excitations are found by decomposing the boson 
field operator <f)j(z,t) around the ground state Hartree 
wave function as 



4>j(z,t)= <j>j(z) + 8<f>j(z,t) 

where the perturbation operators 5<j) 
son commutation relations: 



(6) 



j (z, t) satisfy the bo- 



5cf>j(z,t), 5<f>l(z',t) = 5 jk 5(z 



(7) 



Substituting (||) into (||) and keeping only the leading- 
order terms leads to a time-independent Gross-Pitaevskii 
equation from which the ground-state wave function 
4>j(z) and chemical potential htvo can be determined. 

We work in the condensate rest frame and seek plane 
wave solutions corresponding to a uniform beam of atoms 
moving along the wave guide. The densities of the two 
components are then given by 



Pa = |0a' 2 
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Pb = \<j>b\ = 



Va + Vb- 2V X 



(8) 



where we have assumed a fixed total density p = p a + Pb- 
The relative phase between the two components is ar- 
bitrary, as it is either random, or fixed by the linear 
coupling whose strength we have assumed to be negli- 
gible. The frequency of phase rotation for the ground 
state (chemical potential divided by K) is determined to 
be 



2{V a Vb-V x 2 )p+{Vb-Va)S 
2(V a + V b - 2V X ) 



(9) 



In the case V a = V b , Eqs. (|) and (g) reduce to Eqs. (8) 
and (12) of Ref. respectively. 

The phase-rotation frequency loq of the condensate 
ground state contains two contributions, one being pro- 
portional to the total density p, and the other to the 
detuning S. Both terms result from the combined effects 
of cross- and self-modulation. In contrast to the case of 
a single-component condensate, Eq. (0) shows that the 
cross and self-phase modulation contributions can con- 
spire to cancel the density-dependent term, provided only 
that the condition 



V x 2 = V a V b and V a + V b -2V X ^Q 



(10) 



is met. When the solutions given by Eq. (g) give nega- 
tive densities no homogeneous ground state exists. Thus 
we must verify whether or not these conditions can be 
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fulfilled for positive densities p a and pb- We consider the 
most common situation where all two-body interactions 
are repulsive, V a , Vb, V x > 0, and assume without loss of 
generality that Vb > V x , which implies from Eq. (|l^) that 
V a < V x . It is then easily shown that the requirement of 
p a > and p > yields 



T4 

v b " (v b - V x )p 
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<1, 



(11) 



which can be achieved by an appropriate choice of the 
detuning 8 and/or of the condensate density p. 

To show that these conditions can be met in a real 
system by modifying the trapping potential, we consider 
the two-component 87 Rb condensate, where the compo- 
nents a and b correspond to the hyperfine Zeeman sub- 
levels \F = 2,mf = 1) and \F — 1,TO/ = —1), respec- 
tively. The phase modulation constants are in the ratio 

9a • Qx 

: gb = 0.97 : 1.00 : 1.03. We consider the case of 
harmonic trapping potentials U a = mw (x 2 +y 2 )/2 and 
Ub = muj 2 [(x — x ) 2 + y 2 ]/2, the offset x in their centers 
xo being a control variable that can be changed via a bias 
magnetic field. We have then 



V a b = 
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9a, b 
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where £ = ^Jfi/muj is the extension of the trap ground 
state. Hence to have condition ([To]) satisfied, one can 
choose xo = 0.03£. 

When the condition ( ^0| ) is satisfied, we obtain the 
remarkable result that the evolution of the condensate 
phase is independent of density. We emphasize that by 
itself, this is of course not sufficient to guarantee mea- 
surements unperturbed by nonlinear phase shifts can be 
carried out. Indeed, the detection of any signal relies on 
some departure of the state of the condensate away from 
its ground state. Assuming that this change is small, 
we can describe it in a linearized approach following the 
Bogoliubov treatment. We proceed by expanding the 
Hamiltonian ([!]) with the help of (0) and (^), keeping 
only quadratic terms in the operators 8cj>j , yielding 



j=a,b 



H~ ]T / dz6%(z) 



2m dz 2 



hVji^i 2 



6(f>j(z) 



j=a,b 



E o V i / dz \4>p4>](z)54>}(z) + H.c. 



spinor 4>j(z) satisfies the time-independent nonlinear 
Schrodinger equation. 

The effective Hamiltonian Eq. ([l3|) can be diagonalized 
via a generalized Bogoliubov transformation. As there is 
no confining potential in the z-direction, we expand the 
operators 8<j)j(z) onto plane waves as 



l 3 {z) = (27T)- 1 / 2 / dke lkz c (k) 



(14) 



As a consequence of momentum conservation, only the 
operators Cj(k), cj(fc), dj(—k), and ct(— k) are coupled. 
In order to diagonalize the Hamiltonian ( |l5| ) we introduce 
the annihilation operators for quasiparticles with well- 
defined momentum k according to 

b^k) = J2 [^{k)cj(k) + v M {k)c]{-k)] . (15) 



Invariance under rotation of the coordinate axes clearly 
requires that u^j{—k) = u^j(k) and v^j(—k) = v^(k). 
The coefficients u fjL j{k) and are determined by the 

requirements that the operators & M (fc) and bt(k) obey the 
bosonic commutation relations 



M*)>M*') 
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and 



6 M (fc),St(fc')] = 8 llv 8{k-k'), 
and that the Hamiltonian (]l^) takes the form 

H = J2f dkhcj^k)bl(k)b^k), 



(16) 



(17) 



(18) 



where the oj^k) are thus the frequencies of the elemen- 
tary modes for small collective excitations of the conden- 
sate. 

For a fixed momentum k the coefficients u ll j{k) and 
v ll j{k) may be represented as the matrix elements of the 
2x2 matrices U and V. In order to satisfy the commu- 
tation relations ( |l6|) and ( |l7| ) it is sufficient that 



UU T 



VV f = I, 



UV T = VU T . 



(19) 



The functions u U j(k) and v^j{k) can then be determined 
by substituting Eq.(|15|) into the commutator [b^k), H] = 
Hu) fi(k)b fj,(k) , which guarantees that the Hamiltonian is 
of the form (pi). This yields: 



(20) 



+ W x / dz 0* a (j)l8cj) a (z)54>b(z) + (j) a (f>l8$l(z)84> b (z) where = {u m , u^b, -v^ a , -v^ b ) T and the matrix M 



is given by 



-H.c. 



(13) 



No first-order contributions in 5<f)(z) appear in Eq. ( |l3| ) 
as a consequence of the fact that the condensate 
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where Hj = hk 2 /(2m) + Vjpj. The eigenfrequencies 
u>n(k) whose corresponding eigenvectors satisfy Eqs. (|l9| ) 
e 

w M (fc) 



are 



2m 



hk 2 



2m 



+ V a p a + V b p b 



± [{VaPa ~ Vbpbf + 4V 2 p a p b 



1/2 



1/2 



(21) 



These solutions therefore yield the two branches of the 
quasiparticle excitation spectrum for the two-component 
condensate. 

It is straightforward to show that one of the branches 
does not depend on the density p when the condition (JTT 
is satisfied. In this case the eigenfrequencies are 



w_(Ar) 



Hk 2 

— + 2(V a Pa + VbPb) 

2m 



(22) 



hence w_ corresponds to the 'collisionlcss' branch. 

The existence of a collisionless branch is related to the 
invariance of the system under translation along the z- 
axis. It is easy to show that the corresponding eigenvec- 
tor in this branch is given by 

«7_(fc) = (V x 4> b , ~V a 4>a, 0, Of 

which is fc-independent. The condensate wave function 
is a plane wave, but the choice of inertial frame in which 
this plane wave is at rest (k = 0) is arbitrary. The can- 
cellation of self- and cross-phase modulation is achieved 
by having a uniform density along the z-axis, with a 
constant ratio of the components a and b. If a frac- 
tion of the total atoms are boosted into a moving frame 
while still maintaining their internal superposition state, 
then the relative densities of the two components are not 
changed. This leads to the possibility of imparting kinetic 
energy onto the system without affecting the balance re- 
quired to eliminate density-dependent phase shifts. As 
the coefficients v-j are found to be zero, it is clear that 
the quasiparticles created by b^_(k) correspond simply 
to condensate atoms boosted into a different momentum 
eigenstate. This remarkable property, reflected in the 
collisionless branch of ( |2^ ) should be particularly useful 
for atom interferometry, where the ability to use Bragg- 
pulscs d^j2^] to induce transitions between the conden- 
sate state and the collisionless branch could lead to the 
design of beam-splitters which maintain mean-field-frcc 
conditions. 

In conclusion, we have shown that for certain values 
of collisional constants, the nonlinear phase shift in a 
two-component Bose condensate can be completely elim- 
inated. Such an effect results from the interplay between 
self- and cross-phase modulation between the two com- 
ponents. While we have explicitly shown how it is pos- 
sible to adjust the values of these collisional constants in 



a one-dimensional atomic waveguide, the same method 
also applies to two-dimensional systems with strong con- 
finement in the third dimension. We expect that these 
results will play an important role in atom interferometry 
where uncontrollable nonlinear phase shift limits appli- 
cations in precision measurement. 

Finally, the question remains to determine whether the 
elimination of the mean-field shifts can be achieved in a 
more typical 3-d trap as well. In that case, the ground 
state must be determined from a full three-dimensional 
time-independent set of Gross-Pitaevskii equations. It 
will be interesting to see if the chemical potential can 
be made density-independent in this case by appropri- 
ately engineering the trapping potentials of the different 
components. Because the existence of the collisionless 
branch of the excitation spectrum appears to be related 
to translational invariance, future studies will be required 
to determine whether such a branch exists in the case of 
a 3-dimensional trap as well. 
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